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ABSTRACT 

We detail an innovative new technique for measuring the 2-D velocity moments 
(rotation velocity, velocity dispersion and Gauss-Hermite coefficients I13 and I14) of 
the stellar populations of galaxy halos using spectra from Keck DEIMOS multi-object 
spectroscopic observations. The data are used to reconstruct 2-D rotation velocity 
maps. 

Here we present data for five nearby early-type galaxies to ^3 effective radii. 
We provide significant insights into the global kinematic structure of these galaxies, 
and challenge the accepted morphological classification in several cases. We show that 
between 1-3 effective radii the velocity dispersion declines very slowly, if at all, in 
all five galaxies. For the two galaxies with velocity dispersion profiles available from 
planetary nebulae data we find very good agreement with our stellar profiles. We find 
a variety of rotation profiles beyond 1 effective radius, i.e rotation speed remaining 
constant, decreasing and increasing with radius. These results are of particular impor- 
tance to studies which attempt to classify galaxies by their kinematic structure within 
one effective radius, such as the recent definition of fast- and slow- rotator classes 
by the SAURON project. Our data suggests that the rotator class may change when 
larger galacto-centric radii are probed. This has important implications for dynamical 
modeling of early-type galaxies. The data from this study are available on-line. 

Key words: galaxies: general, galaxies: individual; (NGC 821; NGC 1400; NGC 1407; 
NGC 2768; NGC 4494), galaxies: kinematics and dynamics, galaxies; structure, stellar 
dynamics, techniques; spectroscopic, 



1 INTRODUCTION 

Since the pioneering work of Simkin (1974), major advances 
in the measurement of stellar orbits in early-type galaxies 
have been made. The first quantities to be measured were 
the line-of-sight rotation velocity (V ro t; which can be con- 
sidered a measure of ordered motion within a galaxy) and 
velocity dispersion (a; a measure of the random motions). 
By combining these two measures with photometric ellip- 
ticity (e) it was possible to probe the anisotropy of stellar 
orbits in galaxies (Binney 1978). 

Further progress was made in this field when, in the 
1980s, Bender (1988) showed that although the isophotal 



shapes of early-type galaxies were well described by ellipses, 
small variations from ellipticity could be used to categorise 
them as either 'boxy' or 'discy', and that these variations re- 
late directly to the stellar orbital isotropy. The combination 
of surface photometry and detailed kinematics were there- 
fore demonstrated to provide significant insights into the 
global structure of galaxies. 

More recently, higher order moments of the line-of-sight 
velocity distribution (LOSVD) have been measured (e.g. 
Bender 1990; Rix & White 1992; Gerhard 1993; van der 
Marel & Franx 1993). These characterise deviations from 
a pure Gaussian shape, e.g. skewed and symmetric devia- 
tions, denoted I13 and I14 respectively. Such parameters can 
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indicate kinematic substructures, such as embedded stellar 
discs. 

The development of integral field units (IFUs), such as 
the SAURON instrument, have given new impetus to this 
field, as these provide 2-D maps of the LOSVD. This al- 
lows the definition and estimation of a number of important 
new parameters that can be averaged over the surface of a 
galaxy in a self-consistent way. Perhaps the most important 
of these is the 'spin' parameter (Ah) which is a proxy for 
the projected angular momentum per unit mass. Applied 
to the SAURON galaxy sample of 48 early-type galaxies, it 
was used to show that a dichotomy exists between 'slow-' 
and 'fast-' rotators (Emsellem et al. 2007). Trends between 
photometric and kinematic axis ratios and position angles 
have also been shown to vary with \r; with slow-rotators 
often exhibiting strong misalignments. These observations 
provide new clues to the formation mechanisms of these dif- 
fering galaxy types (Krajnovic et al. 2008; Cappellari et al. 
2007; Thomas et al. 2007). 

Statler & Smecker-Hane (1999) used a different ap- 
proach to building 2-D maps of the LOSVD in NGC 3379. 
Using long-slit spectroscopy along four position angles, they 
were able to reconstruct a model of the 2-D kinematics by 
making suitable interpolations between their observed axes. 
However, such studies are rare due to the large amount of 
telescope time required to observe multiple position angles 
out to large radii. 

Despite these advances, most observations typically 
reach to less than one effective radius, with only a handful of 
studies extending this coverage up to 2-3 R e ff (but see Saglia 
et al. 1993; Graham et al. 1998; Mehlert et al. 2000; Proctor 
et al. 2005; Sanchez-Blazquez et al. 2006; Forestell & Geb- 
hardt 2008). However, such outer regions potentially provide 
key information regarding the physical properties of galax- 
ies. For example, radii within 1 R e fr only sample 50% of the 
baryonic mass and ~ 10% of the angular momentum of an 
idealised galaxy with an R 1 / 4 like light profile and a flat ro- 
tation profile. The situation is even worse when considering 
the distribution of total galaxy mass (including dark mat- 
ter). For instance, the simulations of Dekel et al. (2005) show 
that only ~5% of the total mass is contained within 1 R e fi 
(see also Mamon & Lokas 2005). Furthermore, the signa- 
tures of some formation mechanisms only become apparent 
at large galacto-centric radii (e.g. McMillan, Athanassoula 
& Dehnen 2007; Weil & Hernquist 1996) . A means to probe 
the kinematics in the outer regions of early-type galaxies, 
and hence include a significant fraction of a galaxies' angu- 
lar momentum and mass, is therefore highly desirable. 

In this paper we measure stellar LOSVDs from deep 
multi-object spectroscopic observations and develop a new 
technique for their analysis. Specifically , we measure ro- 
tation velocity, velocity dispersion and the Gauss-Hermite 
coefficients I13 and h.4, out to ~3 R c fr in five nearby early- 
type galaxies. By extending the region studied to such radii 
we encompass more than 85% of the baryonic mass, 30% of 
the angular momentum and 15% of the total mass, a sig- 
nificant increase in the fractions so surveyed in galaxies. To 
achieve this we exploit the large aperture of the Keck tele- 
scope, the high sensitivity of the DEIMOS instrument at 
8500 A and the strong calcium triplet absorption features 
at these wavelengths. Even with this configuration, expo- 




Figure 1. DSS images of the 5 target galaxies. Small red circles 
mark the location of slits from which galaxy background spectra 
were obtained. Large red circles represent 1, 2, 3 and 4 effective 
radii. 



sures of ~ 2 hours are required for our analysis. 

The paper is organised as follows: In Section 2 we out- 
line the observations and data reductions, detailing the tech- 
nique by which we extract the host galaxy spectra. In Sec- 
tion 3 we outline our new kinematic analysis technique. Re- 
sults are presented in Section 4. We discuss our findings in 
Section 5, and present our conclusions in Section 6. 



2 OBSERVATIONS AND DATA REDUCTION 
2.1 Sample selection 

The galaxies in our sample were selected with the intention 
of obtaining spectra of their globular clusters (GCs). The 
galaxies are all nearby, early-types with moderately rich GC 
systems. The selection also had the requirement that reces- 
sion velocities were such that the calcium triplet lines would 
not generally be adversely affected by strong sky lines. The 
main science goal was to obtain GC kinematics, and hence 
probe the mass profile of their host galaxy halos. Initial re- 
sults for NGC 1407 from the 2006 Nov. run have recently 
been reported in Romanowsky et al. (2008). 



© 0000 RAS, MNRAS 000, 1-19 



2-D Kinematics to 3 R e g 3 



2.2 Observations 

The data used in this work are from three observing runs 
using DEIMOS on the Keck-2 telescope. The first observ- 
ing run was carried out in 2006 Nov. using the f 200 1/mm 
grating centered at ~7500 A and a f arcsec slit. This allows 
coverage of the calcium triplet at 8498, 8542 and 8662 A 
at a resolution of ~1 .5 A. The second run was carried out 
in 2007 Nov. and the third in 2008 April, both using the 
same instrument setup as the first run. Seeing was ~0.5 to 

0. 8 arcsec on all nights. A total of five galaxies were ob- 
served (NGC 821, NGC 1400, NGC 1407, NGC 2768 and 
NGC 4494). 

Between two and six independent masks were observed 
per galaxy and the 16 x 5 arcmin field-of-view allowed the 
observation of ^80-120 spectra in each mask. However, only 
~20% of these were at radii from which the signal of the 
background host galaxy was sufficient for the extraction of 
spectra (see Table 4). 

Fig. 1 shows the positions of slits from which host galaxy 
spectra were extracted. We note that it was only possible to 
measure recession velocities in some of the outermost slits, 

1. e. velocity dispersion and higher order moments are not 
presented for some of the outermost positions. 

2.3 Data reduction 

Data reduction were carried out using the DEIMOS DEEP2 
pipeline software. This package produces a variety of out- 
puts from the raw DEIMOS data, including both object and 
'sky' spectra. The data reported here probe the kinemat- 
ics of the stellar populations in the galaxy halos. These are 
based on the 'sky' spectra that were subtracted from the 
globular cluster spectra (the prime targets of the observa- 
tion, see Romanowsky et al. 2008). The 'sky spectra actually 
consist of a co-addition of the true sky and the background 
galaxy light. In order to distinguish them from the true sky 
spectrum we shall therefore henceforth refer to these as the 
'background' spectra. 

The recovery of the galaxy spectra requires an estimate 
of the true sky spectrum in each of the background spectra. 
For each mask, the sky spectra were estimated from the av- 
erage of four to six of the background spectra sampling large 
distances from the galaxy centre (see below). These gener- 
ally lay beyond ~6-7 R c ff . Even at this distance some galaxy 
light is still present. However, since the data reported here 
are limited to within ~3 R e ff, the galaxy light contamina- 
tion of the sky estimates constitutes only a small fraction of 
the galaxy light in our science spectra. Assuming a de Vau- 
couleur's (1953) R 1/4 profile, we estimate that, in points at 
3 Reff, the sky estimate contains only <10% of the galaxy 
light in the science spectra. This falls to ~4% at 2 R e ff . Er- 
rors caused by the residual galaxy light in our sky estimates 
are therefore likely to be much smaller than the quoted er- 
rors on data points at these radii, which are often as much 
as ~30% of the measured parameter values. 

In order to select and characterise the sky spectra, a sky 
index was defined (see Table 2 and Fig. 2) . The central band 
of this index was defined to cover the strong sky emission 
lines between 8600 and 8700 A. The sidebands of the index 
were chosen to cover regions of the background spectra that 
were relatively free of both significant sky emission lines and 
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Figure 2. Definition of the sky index. Upper plot shows wave- 
length regions of the two sidebands and the central band in a 
spectrum with significant galaxy halo signal. Lower plot shows 
the region used in the calculation of the sky index in a scaled sky 
spectrum. The index is defined as a ratio, i.e. excess flux in the 
central band (counts in shown in red) divided by continuum flux 
(counts shown in grey). 

galaxy absorption features. For each background spectrum 
the sky index was then calculated as follows: First, the con- 
tinuum level in the region of the central band was estimated 
by integrating under the linear interpolation of the average 
fluxes in the sidebands. Then the excess flux in the central 
band was calculated as the total flux in the region minus 
the integrated continuum flux. Finally, the sky index was 
calculated as the ratio of the excess central band flux to its 
continuum flux (see Fig. 2; bottom). The index has the prop- 
erty that, as one moves away from the galaxy centre (and 
the galaxy light contribution to the background spectrum 
decreases), the index tends to a constant value (see Fig. 3). 
After a careful visual inspection of the spectra, the back- 
ground spectra to be used as sky estimates were selected 
from this constant region. These were co-added to provide 
a higher signal-to-noise estimate of the sky spectrum. 

The index was used to select candidate sky spectra, 
rather than simply choosing the slits at the greatest dis- 
tance from the galaxy centre, as the index allows identifi- 
cation of spectra suffering from reduction problems not de- 
tected by visual inspection. For the most part, these proved 
to be either slits located near the edges of the mask where 
vignetting effects are large, or slits in which residual globu- 
lar cluster light was still present. Fig. 3 illustrates how the 
sky index allows identification (and elimination) of slits with 
such problems. 

Sky subtraction was carried out by consideration of the 
excess flux in the central band as defined above. First, for 
each background spectrum, we scaled the sky estimate such 
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Galaxy Hubble Distance R e ff Mj( PA p h ot Axis ratio V S y S 

Type (Mpc) (arcsec) (mag) (°) (K band) (km s — *) 



NGC 821 


E6? 


22.4 


50 


-23.6 


30 


0.62 


1729 


NGC 1400 


SA0- 


25.7 


29 


-23.8 


35 


0.90 


574 


NGC 1407 


E0 


25.7 


70 


-24.9 


60 


0.95 


1782 


NGC 2768 


SOi/a 


20.8 


64 


-24.6 


93 


0.46 


1327 


NGC 4494 


El-2 


15.8 


49 


-24.8 


173 


0.87 


1335 



Table 1. Properties of galaxies featured in this work. Distances are from surface brightness fluctuations by Tonry et al. (2001), with the 
distance moduli modified by -0.16 according to Jensen et al. (2003). The values for NGC 1400 and NGC 1407 are taken as the average 
of their individual values. Hubble types are from NED and effective radii are from RC3. Absolute K band magnitudes are calculated 
from 2MASS apparent magnitudes and the distances given in column 3. Photometric position angles (PA p i lot ) and K band photometric 
axis ratios are from 2MASS. Systemic velocities are the values derived from our fitting of central data points (see Section 4). 



Band 



Wavelengths 



Sideband 1 8526.0 - 8536.0 A 
Central band 8605.0 - 8695.5 A 
Sideband 2 8813.0 - 8822.0 A 

Table 2. Band definitions for the sky index used in the selection 
and subtraction of sky spectra. 
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Figure 3. Sky index with radius for one set of NGC 4494 ob- 
servations. Open red circles represent points whose spectra were 
co-added to generate the sky spectrum. Two points clearly exhibit 
aberrant behavior. These are not included in the final sample. 



that the resultant spectrum had the same excess counts as 
the target spectrum. Subtraction of the scaled sky spectrum 
therefore resulted in a spectrum with no excess counts in 
the sky band region, i.e. with counts in the sky band region 
equal to the continuum level in the side bands. This method 
provides a sky subtraction accurate to approximately one 
count over the sky band region. We show examples of sky 
subtracted spectra in Fig. 4 in which it can be seen that 
although sky residuals are present, these are not significantly 
larger than the Poisson noise associated with the strong sky- 
lines. Assuming a dark-sky brightness at Mauna Kea of I ~ 
20 mag. arcsec -2 , then our galaxy halo spectra can reach as 
low as 1-10% of the sky background. We note that this sky 
subtraction technique is similar to that used in Proctor et 
al. (2008) and Norris et al. (2008) in the analysis of optical 
spectra from GEMINI/GMOS. 




3600 8650 8700 
Wavelength (A) 

Figure 4. Sky subtracted spectra of NGC 2768. Red lines show 
the best-fit template derived by the pPXF code during measure- 
ment of velocity moments. Spectra are in order of decreasing 
signal-to-noise broadly covering the range within our data (from 
~40 in the top spectrum down to 7 in the bottom spectrum). 



3 ANALYSIS METHOD 

The analysis presented here involves measurement of veloc- 
ity moments (V, a and the Gauss-Hermite coefficients I13 
and h.4) using the pPXF code of Cappellari & Emsellem 
(2004) - an extension of the pixel fitting code of van der 
Marel (1994). 



3.1 Measurement of kinematics 

The pPXF software works in pixel space, finding the combi- 
nation of template stars which, when convolved with an ap- 
propriate line-of-sight velocity distribution, best reproduces 
the galaxy spectrum. Thirteen template stars were used, 
ranging in spectral type from F8 to Ml. This method of kine- 
matic analysis has the advantage that template-mismatch 
issues are all but eliminated. Operating in pixel space, the 
code also permits the masking out of regions of the galaxy 
spectrum during the measurement. This is of particular im- 
portance in our data as, despite the accurate sky subtrac- 
tion, the spectra still suffer from skyline residuals (see Fig. 
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4). The code also allows suppression of the higher order mo- 
ments I13 and I14 (i.e. fitting only Vobs and a). This was 
found to be useful in order to stabilise the fits in low signal- 
to-noise spectra. 

Errors were estimated from a repeat observation on dif- 
ferent nights of one of the masks of NGC 2768. By binning 
the 60 repeat measurements by observed signal (i.e. counts), 
the rms variations of each of the velocity moments in the re- 
peat measurements were calculated as a function of signal- 
to-noise. The trend of rms with signal-to-noise so obtained 
was then used to make estimates of errors in this, and other 
galaxies. 

As a test of the accuracy of our measurements we com- 
pared them to those of SAURON (Emsellem et al. 2004) 
for the two galaxies in common between the studies. This 
was possible as our data samples 10 precise locations (9 in 
NGC 2768 and one in NGC 821) that are also sampled by 
SAURON. The results of the comparison for a and h3 are 
shown in Fig. 5. A small systematic offset of -12.9 km s _1 
is evident in our velocity dispersion data compared to the 
SAURON values, with a scatter of 7.3 km s _1 . The mea- 
sured I13 values are consistent within their errors. We do 
not compare I14, as errors are generally large with respect 
to our measured values. We leave comparison of Vobs to the 
sections presenting our results. 

During the analysis significant discrepancies were de- 
tected between our results for the a values of NGC 1400 and 
NGC 1407 and the long-slit data of Spolaor et al. (2008a). 
As a consequence, the data of Spolaor et al. (2008a) were re- 
processed through the pPXF software using a large number 
of template stars, resulting in significantly improved mea- 
surements (see Appendix A) . 

The measured velocity moments of all five galaxies are 
available on-line. An example of the available data are pre- 
sented in Table 3. 



3.2 Kinemetry 

Our analysis of these data proceeds by consideration of the 
velocity moments in the space of galacto-centric radius and 
position angle (PA; defined as the angle between the point 
in the galaxy and the galaxy centre from North through 
East) . For the analysis, the data were binned in radial ranges 
(i.e. circular shells) for comparison to the expectations for 
rotating isotropic ellipsoids (or, equivalently, inclined discs). 
Fits of the observed velocity (V b s ) data were carried out by 
least-squares minimisation of the data to the equation: 

Vobs = Vsys + Vrot -COS (4>), (1) 

where V sys is the systemic velocity of the system, V ro t 
is the major-axis rotation velocity and (j> is defined by; 

tan{PA- PA kin ) 
tan(4>) = i - '-, (2) 

with PAkin the position angle of the major rotation axis 
and q the rotation axis-ratio. The values of V ro t , PAkin and 
q were all allowed to vary freely. Values of V sys were taken as 
the value derived from the high-signal-to-noise data of our 
innermost points (i.e. generally the points within 1 R e ff). 
This value was then assumed for all radial bins. 




150 _ 200 

cr (km s 1 ) (this work) 




-0.1 0.1 

h 3 (this work) 



Figure 5. Comparison of velocity dispersion measures (top) and 
h 3 (bottom) with SAURON data (Emsellem et al. 2004). Filled 
symbols arc from NGC 2768, and the open symbol from NGC 821. 
In the plot of cr, the one-to-one (solid) and —12.9 km s — 1 offset 
(dashed) lines are shown. In the plot of I13 the one-to-one line 
(solid) and typical errors (dashed) of ±0.04 are shown. 

In the analysis of the rotation properties of our sample 
galaxies we used two types of radial binning. The first was to 
break the data into 2 or 3 independent radial slices (i.e. into 
circular annuli) , the other was to carry out rolling fits again 
with data binned in circular annuli. Results are presented 
as the major axis values with radii expressed in terms of the 
effective radius (Reff) 1 . Errors on these fits were estimated as 
the rms scatter in the results of 50 Monte Carlo simulations 
of the data using the best-fit parameter values and scatter. 

We were not however able to carry out such an analysis 
on the velocity dispersion data, mainly because we detect lit- 
tle to no variation in this parameter over the radii probed by 
our data. We are therefore not able to measure the axis ratio 
of this velocity moment. However, for two galaxies (NGC 821 
and NGC 2768), the SAURON data of Emsellem et al. 2004 
permit estimation of axis ratios of a in the inner regions. For 
NGC 821 the SAURON data appear to show an axis ratio 
close to 1.0 (despite the ellipticity of the galaxy isophotes). 
On the other hand, in NGC 2768 a appears to have a simi- 

1 Note that R e ff is expressed as the radius of a circle with the 
same area as the ellipse encompassing 50% of the galaxy light. 
Therefore, in galaxies with non-zero ellipticity, a point on the 
major axis at 1 R e g does not lie on, but rather within, the half- 
light ellipse. 
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R.A. Dec. V ba V bs error a a error 113 I13 error I14 I14 error 



NGC 2768 



9:11:31.97 


+60:1:38.7 


1300 


26 


152 


24 


0.14 


0.04 


0.10 


0.05 


9:11:28.47 


+60:2:51.5 


1251 


27 


177 


25 


0.05 


0.04 


0.03 


0.05 


9:11:35.29 


+60:2:03.8 


1284 


5 


192 


4 


0.02 


0.04 


0.02 


0.05 


9:11:38.70 


+60:2:57.4 


1346 


26 


154 


24 


-0.02 


0.04 


-0.01 


0.05 



Table 3. Example of the velocity moment data form this work. Columns are; R.A. and Dec (J2000), observed recession velocity (V Q b s 
in km s — 1 ), velocity dispersion {a in km s _1 ) and Gauss-Hermite coefficients I13 and I14. The full table is available on-line. 
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Figure 6. NGC 2768. Velocity moments are plotted against 
position angle. Data points are coloured according to galacto- 
centric radius (black; R<1 R G ff, red; 1 R c ff <R<2 R G g and grey; 
R>2 R e ff). In the plot of V b s best-fit elliptical models at each 
radius are shown as appropriately coloured lines. The systemic 
velocity (V S ys) is shown as a dashed line. The major kinematic 
axis (PAkjn) is at ~90°. In the plot of a average values are also 
shown as appropriately coloured lines. The values adopted for 
these lines are presented in Table 4. 



lar axis ratio to the isophotes. We therefore adopt these axis 
ratios for the analysis of the velocity dispersion data from 
these galaxies. For the remaining three galaxies (all of which 
have axis ratios >0.8 we assume an axis ratio of 1.0. 

We also experimented with the co-adding of spectra at 
similar effective radii for a measurement. However, this pro- 
cess requires the spectra to be de-redshifted prior to co- 
addition, which complicates the handling of the sky-line 
residuals in the absorption features. In addition the results 
proved almost identical to the values obtained by averag- 
ing individual measurements. On the basis that simplest is 
best, we therefore present a values as individual and radially 
averaged values. 

Our analysis is similar in intent to the 'kinemetry' of 
the SAURON paper by Krajnovic et al. (2006). However, 
our data possess an extremely low sampling density (i.e we 
have on average 60 points inside an area of 24 R^ ff compared 
to SAURON's ~3000 points inside 0.4 R^ ff - a data density 
contrast of 3000:1). Our data are also generally of lower 
signal-to-noise than the SAURON data. Since, as noted in 
Krajnovic et al. (2006), such limitations make it impossible 
to perform a full harmonic analysis, we have restricted our 
analysis to only the low order harmonics and fitted our data 
on circles (or more precisely on annuli) rather than ellipses. 
Unable to perform the full analysis, we have also refrained 
from adopting the full notation used by Krajnovic et al. 
(2006). However we note that the rotation velocity (V ro t) we 
report is most directly related to their ki parameter, while 
the kinematic axis ratio is most directly related to their qkin 
parameter. PAkin has the same meaning in both studies. 



4 RESULTS 

In this section we present the results of our velocity moment 
measurements (V b s , 0, h3 and hi). Errors on individual 
data points were derived from the repeat measurement of 
one mask in NGC 2768 (see Section 3) . Errors on rolling fits 
were estimated by Monte Carlo realisations of the best-fit 
parameters. Note that, for the lowest signal-to-noise spec- 
tra, the higher order moments were either suppressed alto- 
gether or are not presented. Tests on the low signal to noise 
spectra show that changes in recession velocity and veloc- 
ity dispersion caused by the suppression of the I13 and I14 
parameters are ~20% of the quoted errors. We also present 
re-constructions of 2-D rotation velocity maps in this sec- 
tion. 

The results of the fits to independent radial slices 
are presented in Table 4. In NGC 821 and NGC 1400 
two radial ranges were defined (NGC 821: R<1.0 R G ff and 
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Radius (kpc) 
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Radius (R eff ) 

Figure 7. NGC 2768. Radial profiles of the major axis rotation velocity (V ro t), c, kinematic and photometric PA (PA^m and PA prio t) 
and axis ratio (q and 1-e). Note that ellipticity of a is assumed to be 0.54, i.e. the photometric ellipticity. The physical radial distances 
(based on data in Table 1) are given in kiloparsec on the top axis. The a data points (black) are the individual measured values, with 
the line representing a rolling average of 5 points. For all other parameters the black data points represent the values of our fits to 
independent regions (Table 4). The results of rolling fits to the data and their 1-cr errors are indicated as black lines. Green data points 
and lines at small radii are from our analysis of the SAURON data of Emsellem et al. 2004. The red squares in PA and axis ratio are 
the 2MASS photometric values at the radius of the K band 20th magnitude isophote, while the red lines show the R band photometric 
data of Peletier et al. (1990). Overall the agreement with the literature values is very good. 



Galaxy 


No 


<R> 


<R/R e ff> 


<cr> 


Vrot 


axis ratio 


PA kin 






(arcsec) 


(Reff) 


(km s ) 


(km s ) 


(q) 


(°) 


NGC 2768 


21 


31 


0.48 


174±4 


98±3 


0.45±0.06 


90±3 




26 


93 


1.45 


152±4 


178±6 


0.30±0.08 


94±2 




15 


165 


2.57 


123±10 


190±16 


0.35±0.13 


88±6 


NGC 4494 


19 


32 


0.65 


133±3 


58±2 


0.80±0.03 


182±1 




30 


75 


1.53 


111±3 


52±3 


0.50±0.09 


179±3 




33 


143 


2.91 


100±5 


54±5 


0.40±0.20 


170±6 


NGC 1407 


28 


49 


0.70 


261±6 


28±4 


(0.95) 


254±8 




39 


102 


1.46 


240±7 


21±6 


(0.95) 


245±14 




11 


173 


2.48 


198±11 


20±10 


(0.95) 


248±15 


NGC 1400 


7 


31 


1.07 


116±3 


53±3 


(0.9) 


34±4 




11 


68 


2.35 


96±7 


30±7 


(0.9) 


40±14 


NGC 821 


14 


41 


0.82 


149±5 


56±8 


0.30 


33 




17 


78 


1.56 


144±6 


6±8 


(0.30) 


(33) 



Table 4. Kinematic data summary for independent radial slices. The definition of the radial binning is outlined in Section 4. The number 
of slits in each radial bin, their average radius and a are given as well as the results of the rotation velocity fits (Equations 1 & 2). Data 
in brackets are assumed (rather than fitted) values. 
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Figure 8. NGC 2768. Reconstruction of the 2-D rotation velocity map using the rolling fits of Fig. 7. The outer boundary of the 
map has been chosen to reflect the photometric axis ratio of the galaxy (0.46). The colour scale (in km s _1 ) is shown on the right. Also 
shown is the SAURON (Emsellem et al. 2004) velocity map of the central regions of the galaxy. The agreement between the data sets is 
extremely good. 



R>1.0 R cff , NGC 1400: R<1.5 R eff and R>1.5 R eff ). In the 
other three galaxies three ranges were defined (R<1.0 Reff, 
1.0 Reft <R<2.0 R c ff and R>2.0 R cff ). 

We first present our results for NGC 2768, as this galaxy 
possesses high signal-to-noise data, a large sample size, and 
is also one of the galaxies observed by SAURON. 



4.1 NGC 2768 

The brightest galaxy in a loose group, NGC 2768 is classified 
as an E6 galaxy in RC3 (de Vaucouleurs et al. 1991). How- 
ever, in both the Revised Shapley-Ames Catalogue of Bright 
Galaxies (Sandage, Tammann & van den Bergh 1981) and 
the Carnegie Atlas of Galaxies (Sandage & Bedke 1994) it is 
classified as an SO1/2 due to a 'definite outer envelope'. The 
galaxy is highly elliptical with K band photometric axis ra- 
tio of 0.46, and PA of 93° ( Jarrett et al. 2000). It also reveals 
a polar orbiting dust and gas ring, possibly of external ac- 
cretion origin (see Martel et al. 2004 and references therein). 

Photometric profiles out to 200 arcsec are presented in 
Peletier et al. (1990). They find a rapidly rising ellipticity 
profile with radius (i.e. a falling axis ratio), and report indi- 
cations of a 'disk-like component' at about 50 arcsec which 
gives way to boxy isophotes beyond about 100 arcsec. 

Fried & Illingworth (1994) probed the kinematics along 
the minor and major axes out to about 50 arcsec. They 
found evidence for stellar rotation aligned along the photo- 
metric major axis and a 'spheroidal' velocity field. They also 
detected minor axis rotation associated with the polar ring. 

The central regions (~80x40 arcsec) of the galaxy were 
also observed with the SAURON integral field unit (Em- 
sellem et al. 2004). The SAURON results show strong rota- 
tion with a 'cylindrical velocity field' and a central dip in a 
(McDermid et al. 2006). They find the galaxy to be a fast- 
rotator, but, contrary to Sandage & Bedke (1994), conclude 
that the galaxy '...does not have evidence for a disc'. They 



therefore adopt the RC3 classification of this galaxy as an 
E6. 

The SAURON team also investigated the very central 
regions (~ 4x4 arcsec) using the OASIS integral field spec- 
trograph (McDermid et al. 2006), identifying emission as- 
sociated with the polar disk and a young (2.5 Gyr) cen- 
tral stellar population. However, other spectroscopic age de- 
terminations in the literature generally prove inconsistent, 
with studies reporting the central stellar population to be 
both old and young (e.g. Howell 2005; Denicolo et al. 2005; 
Sil'Chenko 2006) . It is unclear whether this reflects inconsis- 
tencies in the determination of the properties of the central 
populations, or is rather a real effect, perhaps caused by 
variable sampling of the central star-forming disk. 

Hakobyan et al. (2008) suggest this galaxy may be a 
merger remnant with the presence of a young stellar pop- 
ulation suggested by "...the existence of strong HI and CO 
emission, and the presence of dust and ionised gas." . The 
recent type lb supernova SN2000ds (Filippenko & Chornock 
2000) supports this conclusion. 

Presenting our results, we first consider the distribution 
of the NGC 2768 velocity moments with PA. Results are 
plotted in Fig. 6 with point colour depicting galacto-centric 
radius. Values derived for each radial bin are presented in 
Table 4. 

The rotation in this galaxy can be seen to be strongly 
elliptical. An increasing rotation velocity with radius is also 
evident. The I13 parameter reflects this behaviour, with near- 
zero values in the slowly rotating inner regions (R<1 R ff). 
However, at larger radii I13 increasingly mirrors the elliptical 
rotation curve. The increasing signal in I13 with radius is also 
detected in the SAURON maps of the galaxy (Emsellem et 
al. 2004). 

Both a and I14, however, exhibit unusual behaviour, as 
they do not possess the minor-axis symmetry expected from 
even moments (see Krajnovic et al. 2008). This must be 
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understood in terms of the degeneracy between a and lu 
discussed at length in Cappellari & Emsellem (2004). The 
degeneracy is worst when; i) signal-to-noise is low, and ii) 
when two distinct kinematic components have roughly equal 
luminosities. To test that this was indeed causing the lack 
of axi-symmetry in our data, we re-fitted the data for this 
galaxy using pPXF with I13 and lu suppressed (i.e. set to 
zero). The test confirmed that, if I14 is forced to zero, the 
dip in a becomes visible on both sides of the galaxy. 

Given that the I14 parameter is more prone to the ef- 
fects of low signal-to-noise and sky-line residuals than a, we 
conclude that, in its outer reaches, this galaxy probably ex- 
hibits a significant drop in a on both sides of the major axis, 
as would be expected from the disc of an SO galaxy. 

Velocity dispersion appears to be decreasing with in- 
creasing radius (Fig. 6 and Table 4). However, the dip on 
the major axis (PA~270°) is clearly effecting measured aver- 
age values, particularly since the slits in the outer regions of 
the galaxy from which we were able to recover a galaxy spec- 
trum tend to lie along the major axis where the galaxy light 
is brightest (Fig. 6). The gradient in a presented here may 
therefore be somewhat exaggerated by this effect. We note 
for future reference that the four data points with a <100 
km s _1 all lie very close to the major rotation velocity axis. 

To measure variations with radius we carried out a 
rolling fit to the recession velocity data in annuli at increas- 
ing radii allowing V ro t, PAkin and kinematic axis ratio (q 
in Equation 2) to vary with radius. The fits were performed 
using a minimum of 10 data points, increasing to 20 in the 
outer regions. Results of these fits are presented in Fig. 7. We 
also performed rolling fits to the SAURON data (Emsellem 
et al. 2004) for comparison. 

The agreement in V ro t between our results and the 
SAURON data is good, with near identical rotation veloci- 
ties over the whole ~20 arcsec overlap region. However, the 
rotation velocity clearly continues to rise beyond the extent 
of the SAURON data, with the rotation velocity doubling 
to ~200 km s" 1 between 1 and 3 R ff- This has a signifi- 
cant impact on the value of V/a obtained for this galaxy, 
underlining the importance of data at large radii in the con- 
struction of even simple kinematic descriptions of galaxies 
(see Section 5). 

Also shown in Fig. 7 is a plot of a with radius. In this 
plot the agreement between the SAURON results (Emsellem 
et al. 2004) and ours can be seen to be reasonable. However, 
the comparison of a presented in Fig. 5 is largely based on 
data from this galaxy, and exhibits a ~13 km s~ x systematic 
offset. The results indicate a strong radial trend in a within 
1 R e ff, which levels off somewhat beyond. However, we again 
note the increasing tendency for our outer data points to lie 
along the major kinematic axis and the trend for such points 
to possess depressed a. To get some idea of the impact these 
trends might have on our results, consider the four points 
with the lowest a (<100 km s _1 ), which, as noted in the 
previous section, all lie very close to the major kinematic 
axis. If these points are excluded from our considerations 
we would conclude that a was flat beyond ~1 R e ft- 

We detect a PAkin of 91±3°throughout the galaxy in 
good agreement with the photometric PA from 2MASS (Jar- 
ret et al. 2000) and Peletier et al. (1990). We detect no vari- 
ation in PAkin with radius, with our measurements of the 
SAURON data yielding the same, near-constant value. 
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Figure 9. NGC 4494. Same as Fig. 6. The major axis PA is 
173° 

Comparison of our rotation axis ratio measurements to 
those of SAURON again give good agreement in the over- 
lap region. The results clearly indicate that the rotation in 
the central regions, variously described as 'cylindrical' (Em- 
sellem et al. 2004) or 'spherical' (Fried & Illingworth 1994), 
rapidly give way to much more disk-like values at larger 
radii. Interestingly, the comparison with the photometric 
axis ratios of 2MASS and Peletier et al. (1990) show the 
rotation axis ratio to pass from less than the photometric 
axis ratio in the outer regions to greater than the photo- 
metric axis ratio in the inner regions. The most natural ex- 
planation for the galaxy having a smaller photometric than 
kinematic axis ratio in the inner regions is the presence of a 
highly inclined, but inflated disc, such as would be found in 
an SO galaxy, in accord with the Sandage et al. (1981) and 
Sandage & Bedke (1994) classifications. 

To aid in the visualisation of our results, and to facili- 
tate the comparison of our results with those of SAURON, 
we have reconstructed the 2-D velocity map of NGC 2768 in 
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Figure 10. NGC 4494. Same as Fig. 7 with the following exceptions. Green data points in V ro t and <r plots are the major axis data 
from Napolitano et al. (2009). The planetary nebula data (blue stars) are also from Napolitano ct al. (2009). Grey data points are the 
data of Bender, Saglia &: Gebhardt (1994). Red lines in plots of PA and axis ratio are the V band results of Napolitano et al. (2009). 
Large offsets between the kinematic and photometric axis ratio and position angles can be seen. 



Fig. 8. The image was generated using the rolling fit values 
of Vrot , PAkin and rotation axis ratio of Fig. 7. Agreement 
with the SAURON results can be seen to be extremely good. 
Also evident is the progression from 'cylindrical' rotation in 
the inner, bulge dominated regions to discy rotation in the 
outer regions probed by our data. 

In summary, our results for this galaxy are all in accord 
with the Shapley-Ames Catalogue (Sandage et al. 1981) and 
Carnegie Atlas of Galaxies (Sandage & Bedke 1994) classifi- 
cations of this galaxy as an SO, with the disk-like kinematics 
in the outer regions giving way to the more 'cylindrical' ro- 
tation in the inner regions. 



4.2 NGC 4494 

NGC 4494 is an elliptical galaxy in the Coma I cloud (Forbes 
et al. 1996) with a K band photometric position angle of 
173°and axis ratio of 0.87 (Jarret et al. 2000). It reveals a 
small central dust ring with ~60 degree inclination (pho- 
tometric axis ratio ~ 0.5) aligned with the major axis of 
the galaxy (Carollo et al. 1997). The central regions also 
host a kinematically distinct core (Bender 1988). The mass 
modeling analysis of Kronawitter et al. (2000) indicates a 
rising M/L profile. However, evidence for dark matter in 
NGC 4494 has been questioned by Romanowsky et al. (2003) 
based on an extended a profile from planetary nebulae kine- 
matics (but see also Dekel et al. 2005 and Napolitano et 
al. 2009). The stellar population analysis of Denicolo et al. 
(2005) found an intermediate age of 6.7 Gyr in the central 



regions, indicating that this is not a simple, old, passively 
evolving elliptical galaxy. 

Our analysis of this galaxy was the same as that for 
NGC 2768. The distribution of the velocity moments with 
PA are shown in Fig. 9. Results show the rotation in this 
galaxy to be strongly elliptical, but rotation velocity shows 
no sign of variation with radius. Within 1 R e ff, the I13 pa- 
rameter shows the anti-correlation with V ro t expected from 
a disc embedded in a more slowly rotating (or non-rotating) 
bulge, although it is not possible to trace this behaviour out 
to the outer regions where the signal-to-noise is low. We find 
a decreasing a with radius which flattens out beyond 1 R crr . 
We find no evidence for variation of a with PA. 

The rolling fits to V b s are presented in Fig 10. These 
fits were performed using 10 data points in the inner re- 
gion rising to 20 data points in the outer regions of the 
galaxy. For V ro t, we find good agreement with profiles in 
the literature. Our results are also consistent with the plan- 
etary nebulae data of Napolitano et al. (2009) in the outer 
regions (notwithstanding their large uncertainties). As de- 
duced from Fig. 9, the galaxy exhibits no sign of variation in 
rotation speed with radius over the whole radial range from 
1-3 R e ff. The reconstructed 2-D velocity map of this galaxy 
is shown in Fig. 11. 

Agreements with the a profiles and planetary nebu- 
lae data in the literature are again good, with a declining 
steadily with radius within 1 R orr , but leveling off beyond. 
Indeed, beyond 1 R cr r , the data are consistent with no slope 
in a at the 1-a confidence level. 
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Figure 11. Reconstructed 2-D velocity maps of (top to bottom) 
NGC 4494, NGC 1407 and NGC 1400. Colour scales (in km s" 1 ) 
are shown to the right of each image. Note that the variations in 
Vrot with radius in NGC 1407 that causes the unusual patterning 
in the reconstruction are only detected at ~1-ct significance. 



The kinematic position angle exhibits a shift of ~ 10 °at 
~1.5 R c ff , mirroring the shift in photometric position angle 
(Napolitano et al. 2009). 

The rotation axis ratio is similar to the photometric axis 
ratio in the innermost regions, but falls rapidly at the radius 
at which the position angle changes. The photometric axis 
ratio (Napolitano et al. 2009), on the other hand, remains 
near constant with radius. The disparity in axis ratios is all 
the more curious when the lack of photometric evidence for 
an embedded disc at these radii (Carollo et al. 1997) is con- 
sidered. However, it is interesting to note that the kinematic 
axis ratio in the outer regions of the galaxy is similar to the 
photometric axis ratio of the central embedded disc of ~0.5 
(Carollo et al. 1997). This suggests that the inner disc could 
in fact be part of a much larger (warped) disc structure. 

The interpretation of the kinematics in this galaxy is 



clearly complex. The highly elliptical V ro t and I13 at ~1 R e ff 
and their possible association with the central disc would 
seem to suggest that this galaxy harbours an embedded 
stellar disc at large radii. However, the lack of photometric 
evidence for such a disc raises serious questions about such 
an interpretation. An alternative description is suggested by 
the results of Balcells (1991) who shows that the signature 
of rotation accompanied by an I13 signal can be generated 
by a minor merger event, as a result of the re-distribution 
of the angular momentum of the merging galaxy. Such an 
event could also be responsible for the shift in photomet- 
ric and kinematic position angles at ~1.5 R c ff. Our results 
are therefore at least qualitatively consistent with the idea 
that this galaxy is a merger remnant. The presence of a 
kinematically distinct core (Bender 1988) also supports this 
interpretation. 

4.3 NGC 1407 

NGC 1407 is the brightest group galaxy in a dwarf galaxy 
dominated group (Trentham et al 2006). The galaxy has 
been classified as a weakly rotating E0 (Longo et al. 1994) 
with a "core-like" central luminosity profile (Lauer et al. 
1995; Spolaor et al. 2008a). The stellar population analy- 
sis of Spolaor et al. (2008b) found the galaxy to possess a 
uniformly old age within ~0.6 R c ff. 

Plots of velocity moments with PA are presented in Fig. 
12. The data were not fitted for kinematic axis ratio due 
to low signal-to-noise and the large uncertainties in fits of 
axis ratio when the value approaches 1.0. However, it is clear 
from the plot of Voba with PA that the rotation has very low 
ellipticity (kinematic axis ratio ~1.0), similar to the 2MASS 
photometric axis ratio of 0.95. This value was therefore as- 
sumed for the rolling fits. A slow (~25 km s" 1 ) rotation is 
evident in all three radial bins. There is a suggestion of non- 
zero I13 in the inner regions, consistent with the results of 
the re-analysis of Spolaor et al. (2008a) data (see Appendix 
A). The data suggest a declining a with radius. 

Rolling fits were performed on the V b s data using 12 
data points in the inner region rising to 20 data points in the 
outer regions of the galaxy (Fig. 13) . The results confirm the 
near constant V ro t with radius between ~0.5 and 3 R e e (we 
note that the variations with radius are only of order 1-er 
significance) . Fig. 13 also confirms the shallow decline in a 
with radius identified above. The reconstructed 2-D velocity 
map of this galaxy is shown in Fig. 11. 

PAkin shows a ~2-ct deviation from the photometric PA 
of Goudfrooij et al. (1994) in the inner regions (< 1.0 R c ff), 
but is in good agreement with the 2MASS photometric PA 
at 1.5 Rcff. The cause of the apparent discrepancy in the 
inner regions is unclear, although significant triaxiality is 
one possible cause. However, it may also simply be the re- 
sult of the large uncertainties that arise when fitting both 
photometric and kinematic data when axis ratios approach 
unity. 

4.4 NGC 1400 

Variously classified as an E (e.g. Da Costa et al. 1998) and 
an SO (e.g. 2MASS; Jarrett et al. 2003), this galaxy is asso- 
ciated, but apparently not interacting with, the NGC 1407 
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Figure 13. NGC 1407. Same as Fig. 6 with the following exceptions. Kinematic axis ratio is not presented as the proximity of the 
rotation axis ratio to 1.0 makes its constraint beyond the capabilities of the signal-to- noise of our data. In the Wot and a plots green 
points are from our re-measurement of the data of Spolaor et al. (2008a) (see Appendix A). The red line in the PA plot is the I band 
data of Goudfrooij et al. (1994). Photometric PAs are shown with a 180 degree offset for the purpose of comparison to the kinematic 
data. 



group. Spolaor et al. (2008b) found the stellar population to 
possess a uniformly old age out to ~1.3 R c ft. 

The distribution of V bs and a with PA are shown in 
Fig. 14. As for NGC 1407, the kinematic axis ratio was not 
fitted. Instead, the 2MASS photometric axis ratio of 0.9 was 
assumed. For this galaxy, the signal-to-noise was too low for 
the accurate measurement of the I13 and I14 moments. These 
were were therefore suppressed during kinematic analysis 
using the pPXF code (Section 3). However, we note the 
significant I13 detection in the re-analysis of the Spolaor et 
al. (2008a) data (Appendix A). Our data sampled only one 
data point within 1 R ff. Radial binning was therefore car- 
ried out in the bins R<1.5 R c ft and R>1.5 R e ff. For some of 
the outermost points errors in velocity dispersion exceeded 
100 km s _1 . These data were omitted from both plots and 
analysis. 

Fig. 14 shows that V ro t is lower in the outer radial range 
than in the inner range. There is also the suggestion of a 
falling a with radius. The data were not fitted for rotation 
axis ratio due to low signal-to-noise and the large uncer- 
tainties in fits in axis ratio measurements when the value 
approaches 1.0. However, it is clear from the plot of Vobs 
with PA that the rotation has very low ellipticity (rotation 
axis ratio ~1.0), similar to the 2MASS photometric axis ra- 
tio of 0.90. This value was therefore assumed for the rolling 
fits. 

Due to the small sample size, rolling fits were performed 
to Vobs using 7 points in the inner region rising to 12 in 
the outer regions. The resultant plots with radius (Fig. 15) 
confirm the conclusions drawn from the plot with PA of a 
falling Vrot and a with radius. Indeed, this galaxy shows the 
sharpest fall-off in a with radius in our sample - falling from 



>300 km : 



in the centre to ~100 km s at ~1.5 R c ft, but 



remaining near constant beyond. However, there appears to 
be a significant unexplained offset between a from this work 
and the values from the data of Spolaor et al. (2008a) in the 
overlap region (both before and after the re-measurement 
of the data using the pPXF; See Appendix A). It is impor- 
tant to note that this discrepancy increases with decreas- 
ing radius, with the discrepancy rising from ~ 20 km s 
at 1.4 R c ff to 60 km s" 1 at 1 R e ft (see also Fig. 5). This 
clearly indicates that this is not a signal-to-noise issue. A 
detailed examination of the spectra of the two studies of 
NGC 1400 could identify no obvious reasons for such a dis- 
crepancy. Consequently, we must contemplate the possibility 
that this is a real effect. We therefore note that the sense of 
the discrepancy is consistent with the idea that the optical 
and NIR probe different stellar populations, with the NIR 
being dominated by a lower velocity dispersion (more cen- 
trally concentrated) population than the optical. The strong 
stellar population metallicity gradient in NGC 1400 (-0.38 
dex/dex; Spolaor et al. 2008b) may therefore be the cause of 
the large discrepancy in this galaxy. The reconstructed 2-D 
velocity map of this galaxy is shown in Fig. 11. 

PAkin shows no significant variation with radius and is 
in good agreement with the photometric PA of 2MASS and 
Spitler et al. (2009; in prep). 

With regard to this galaxy's rather uncertain morpho- 
logical classification, the evidence for a disc with a declining 
rotation, but near constant a profile would seem to favour 
its classification as an elliptical with an embedded disc. 

4.5 NGC 821 

NGC 821 is an isolated (Reda et al. 2007), elliptical galaxy 
with morphological classification E6. The galaxy is con- 
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Figure 12. NGC 1407. Same as Fig. 6 except that the rotation 
axis ratio was not fitted but assumed to be equal to the 2MASS 
photometric axis ratio (i.e. 0.95). The major axis PA is 60°. 



sidered to possess a disk embedded in a spheroidal halo 
(Michard & Marchal 1994; Ravindranath et al. 2002; Em- 
sellem et al. 2004). The photometric profiles of Goudfrooij 
et al. (1994) show that ellipticity peaks around 15 arcsec 
(~0.3 Reft). Kinematic profiles have been widely published 
(e.g. Pinkney et al. 2003; Proctor et al. 2005; Forestell & 
Gebhardt 2008), although there is still some debate as to 
whether the a profile is falling or not and whether NGC 821 
contains dark matter (Romanowsky et al. 2003). The stellar 
population study of Proctor et al. (2005) reveals a young 
central stellar population which gives way to old ages by 1 
effective radius. 

Plots of velocity moments with PA are shown in Fig. 
16. Only the data within 1 R c ff were fitted for position angle 
and rotation axis ratio due to low signal-to-noise and small 
sample size in the outer regions. The values derived from 
the inner regions was then assumed for the outer regions. 

Fig. 16 shows strong rotation within 1 R e ff, but exhibits 
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Figure 14. NGC 1400. Same as Fig. 6 with the following ex- 
ceptions. The velocity moments I13 and I14 were suppressed in this 
low signal-to-noise sample. Radial bins of the independent data 
points are R<1.5 R e fj and R>1.5 R e ff- The rotation axis ratio 
was not fitted but assumed to be 0.9 (i.e. equal to the 2MASS 
photometric axis ratio). 




100 200 
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Figure 16. NGC 821. Same as Fig. 14 except that the radial 
binning is in ranges of R<1.0 R c fj and R>1.0 R c ff only. The major 
axis PA is 30° . 



a value consistent with zero rotation beyond this region. We 
also detect no variation in a with radius. 

Rolling fits of the rotation velocity data (Fig. 17) were 
performed using 7 data points in the inner region rising to 
12 data points in the outer regions of the galaxy. Fig. 17 
confirms the falling V ro t with radius and flat a identified in 
the plots with PA. The figure also show that, where they 
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Figure 15. NGC 1400. Same as Fig. 7 with the following exceptions. The independent data points are for R<1.5 R e ff and R>1.5 R c ff . 
Green data points in V ro t and a plots are from the re-measurement of the data of Spolaor et al. (2008a) (see Appendix A). The red line 
in the PA plot is the Sloan i band data of Spitler et al. (2009; in prep). 



overlap, there is good agreement between our results and 
both the (Emsellem et al. 2004) and Forestell & Gebhardt 
(2008) results for V ro t, o", PAkin and kinematic axis ratio, 
although we note one discrepant data point of Forestell & 
Gebhardt at 1.8 R e ff- 

Our a values are also in good agreement with literature 
studies of both stellar and planetary nebulae kinematics. 
We note that Coccato et al. (2009) also measure a planetary 
nebula rotation velocity consistent with zero for a data point 
at 5 Reff with rotation and a of 14±25 km s _1 and 51±18 
km s , respectively. 

The single kinematic axis ratio value that we derive is 
significantly lower than the photometric axis ratio, confirm- 
ing the result found using the SAURON data (Krajnovic 
2008). 

It is also interesting to note that the extent of the disc 
evident in the V ro t plot of Fig. 17 is similar to the extent 
of both the dip in photometric axis ratio of Goudfrooij et 
al. (1994) (bottom plot; Fig. 17) and to the extent of the 
young central population identified in Proctor et al. (2005) , 
both of which are evident out to nearly 1 R c ff. This lends 
credence to the speculation that the young central ages are 
associated with the presence of a disc that has undergone a 
recent burst of star formation. Also of interest is the fact that 
the photometric axis ratio of the galaxy remains fairly low 
(~0.7) even at radii well beyond the extent of the disc. This 
suggests that the galaxy is comprised of a disc embedded in 
an highly elliptical bulge in agreement with the photomet- 
ric classification in the literature (Michard & Marchal 1994; 
Ravindranath et al. 2002; Emsellem et al. 2004). 

To aid in the visualisation of our results, and to facili- 
tate the comparison of our results with those of SAURON 
(Emsellem et al. 2004), we have reconstructed the 2-D rota- 
tion map of the galaxy in Fig. 18. The image was generated 



using the rolling fit values of V ro t of Fig. 17 and assuming 
rotation to have constant PA^m and axis ratio, with values 
equal to the values within 1 R c h (33° and 0.3 respectively; 
Table 4) . The SAURON map is also shown. Despite the spa- 
tial gap between the two data sets, agreement can be seen 
to be good. 



5 DISCUSSION 

We have demonstrated that our analysis gives good agree- 
ment with the SAURON (Emsellem et al. 2004) study and 
other data from the literature in regions where they over- 
lap. We have also shown that the results provide valuable 
insights into individual galaxies. We next consider our re- 
sults in a more general sense. 



5.1 Rotation 

The five galaxies in our sample exhibit a range of rota- 
tion profiles. Two galaxies (NGC 1407 and NGC 4494) ex- 
hibit flat rotation profiles (within errors) beyond 0.2 and 
0.4 RcH respectively. One galaxy (NGC 2768) exhibits rota- 
tion that increases in amplitude over the full extent of our 
data (i.e. to ~3 R ff). The remaining two galaxies (NGC 821 
and NGC 1400) exhibit rotation curves that decline beyond 
1 Reff, and in the case of NGC 821, becomes consistent with 
zero rotation by 1.5 R c ff. It is therefore evident from our 
data that studies confined to the central regions of galaxies 
can miss vital aspects of their rotation properties. 
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Figure 17. NGC 821. Same as Fig. 7 with the following exceptions. Radial binning is in ranges of R<1.0 R e fj and R>1.0 R c ff only. 
Kinematic PA and axis ratio were only estimated for the inner range. Inner range values of these parameters were assumed for the 
measurement of V ro t- Green data points are from SAURON (Emsellem et al. 2004). Grey data points are from Forestell & Gebhardt 
(2008). Blue stars are the planetary nebula data of Coccato et al. (2009). Red lines in PA and axis ratio plots are the photometric data 
of Goudfrooij et al. (1994). 



5.2 Velocity dispersion 

The five galaxies in our sample exhibit only shallow slopes in 
a with radius beyond 1 R c ff. Indeed, in most cases, results 
in these regions are consistent with zero slope at the 2- 
o" significance level. Such flat a and rotation profiles (with 
the possible exception of NGC 821) represent the classic 
signature of a dark matter halo in each galaxy. These data 
will therefore prove invaluable to future studies that perform 
full dynamical modeling of galaxies. 



5.3 Higher order moments 

We note that all 5 galaxies in this study exhibit non-zero I13 
values, the sense of which mirror the rotation velocity. This 
trend is identified in the NIR spectra of this work in 3 of the 
5 galaxies. In NGC 1400 it is identified in the re-analysis of 
the Spolaor et al. (2008a) optical data, while in NGC 821 it 
is evident in the optical data of Forestell & Gebhardt (2008) . 
When I13 follows the inverse trend of the rotation velocity 
it is often taken as the signature of a cold (rotating) stellar 
population embedded within a hotter (pressure supported) 
system. However, Balcells (1991) shows that this signature 
can also be induced by the action of a merger event. We 
therefore conclude that all five of the observed galaxies are 
either merger remnants or harbour embedded discs. 

There do not appear to be any clear trends in the I14 
parameter in our current data, we suspect that this is due 
to the low signal-to-noise of the majority of our data. 



5.4 Kinematic and photometric misalignment 

For fast-rotators, the SAURON survey revealed that the 
kinematic and photometric position angles are usually well 
aligned, and that their axis ratios are similar (Emsellem et 
al. 2007). This is true for two of our galaxies (NGC 2768 and 
NGC 1400). However, in both NGC 4494 and NGC821 there 
appear to be significant offsets between kinematic and pho- 
tometric axis ratios, while for NGC 821 an offset in axis ratio 
is detected This result for NGC 821 confirms the SAURON 
result, despite the fact that we have only one kinematic data 
point at large radii. 



5.5 Anisotropy diagram 

Binney (1978) showed the usefulness of the anisotropy dia- 
gram (V/c vs photometric ellipticity) as applied to elliptical 
galaxies. In Fig. 19 we show this diagram for the two galax- 
ies in common with the SAURON survey (i.e. NGC 821 
and NGC 2768), with the lines indicating the radial vari- 
ation. The plot shows that the value of V/cr in NGC 821 
declines with increasing radius, while NGC 2768 increases. 
Thus NGC 821 trends away from the isotropic rotator line, 
while NGC 2768 trends towards it. It is therefore clear that 
the central regions of galaxies do not necessarily reflect the 
properties of galaxies as a whole, and care must be taken 
when drawing conclusions from central data only. We note 
that similar conclusions were reached from PN kinematics 
by Coccato et al. (2009). 
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Figure 18. NGC 821. Reconstruction of the rotation profile using the rolling fits of V ro t (Fig. 17). However, beyond 1 R e ff the values 
of kinematic PA and axis ratio are assumed to be the same as within 1 R e ff (33°and 0.3 respectively). The outer boundary of the map 
has been chosen to reflect the photometric axis ratio of the galaxy (0.62). The colour scale (in km s _1 ) is shown on the right. Also shown 
is the SAURON (Emsellem et al. 2004) map of the central regions of the galaxy. The agreement between the data sets is again extremely 
good. 



5.6 The SAURON X R parameter 

As part of the analysis of the SAURON 2-D kinematics, 
Emsellem et al. (2007) defined a new parameter, Ah as; 



Xr = 



< R.\V\ > 



(3) 



where <R.|V|> is the luminosity weighted average of 
the product of radius and rotation velocity over the surface 
of the map and < R.^/iV 2 + er 2 )> is a similarly luminosity 
weighted average. Thus Xr is a measure of the projected 
angular momentum per unit mass (i.e. the specific angular 
momentum. 

Emsellem et al. (2007) found that galaxies were clearly 
distinguished into two classes on the basis of this parameter, 
i.e. 'fast-rotators' (Xr > 0.1) and 'slow-rotators' (Xr < 0.1). 
Fast rotators tended to be lower luminosity, discy ellipticals, 
while the massive boxy ellipticals were often slow rotators. 

The two galaxies in common with the SAURON project 
(NGC 821 and NGC 2768) were both classified as 'fast- 
rotators' (i.e. Xr > 0.1) by Emsellem et al. (2007). Their 
Xr values are very similar (0.258 and 0.268 respectively). 
However, we find that when probed to radii greater than 
1 Reff NGC 821 reveals a declining rotation velocity (and 
roughly constant velocity dispersion) profile. This suggests 
that if the outer region data were included the calculated 
Xr value would be smaller. Although due to the luminosity 
weighting in the Xr definition more emphasis is given to the 



central regions, so the actual change in the value might be 
small. 

To test this, we constructed a model of each of these 
galaxies. Each model assumed R 1 ^ 4 de Vaucouleur's profiles. 
For each galaxy, the rotation parameters (V ro t, PAkin and 
kinematic axis ratio) were then characterised by a polyno- 
mial fit to the plots with radius (Figures 7 and 17) . The lumi- 
nosity and velocity maps so produced were then combined as 
per Emsellem et al. (2007) to calculate the X R values inside 
ellipses of increasing effective radius 2 . The results are pre- 
sented in Fig. 20. Agreement with the published SAURON 
values for the inner regions can be seen to be extremely 
good. However, for NGC 821 our model suggests that, de- 
spite its classification as a fast-rotator (when only data from 
within 1 R c ff is considered), the decreasing rotation at large 
radii could indeed result in the galaxy being classified as a 
slow-rotator when all the data out to 3 R e ft are considered. 

In contrast to NGC 821, NGC 2768 has a profile that 
continues to rise all the way out to ~3 R c ff. Its A_r value 
will therefore continue to increase as larger and larger radii 
are included in the calculation, and the galaxy will remain a 
'fast-rotator'. Therefore, despite their similarity within one 
effective radius, both of these galaxies have vastly different 



2 Although precise values of the sersic index for the two galaxies 
are not available, we note that the results are relatively insensitive 



to this parameter 
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Figure 19. Anisotropy diagram. V ro t/o" is plotted against cllip- 
ticity at various radii for NGC 821 (black) and NGC 2768 (red). 
The remaining three galaxies are not shown as their V ro t / c and 
photometric ellipticities show little or no variation with radius. 
The thin lines represent the inner SAURON data, and the thick 
lines are our data out to 3 R e ff • The smooth curved line represents 
an ideal isotropic rotator. 

distributions of angular momentum with radius (see also 
Coccato et al. 2009). 

Although not observed in the SAURON survey, we as- 
sume NGC 1407 to be a slow-rotator, as its V/a value is low 
(we estimate 0.09 on the basis of Equation 23 of Cappellari 
et al. 2007) and its kinematics are similar to NGC 5982 
which was observed by SAURON (Krajnovic et al. 2008) 
and defined as a slow-rotator (Emsellem et al. 2007). fn the 
case of NGC 1400 and NGC 4494 their relatively high rota- 
tion velocities and relatively low velocity dispersions clearly 
indicate them to be fast rotators. 

5.7 The RMS velocity parameter 

In addition to the V rot /(J and Xr parameters which probe 
orbital anisotropy and angular momentum respectively, an- 
other fundamental parameter of a galaxy is its internal kine- 
matic energy. The observational measure of the total kinetic 
energy is the RMS velocity parameter, which locally is equal 
to the quadrature combination of rotation and dispersion ve- 
locities, i.e. U rms = \/V r 2 ot + a 2 (e.g. Napolitano et al. 2009). 

In Fig. 21 we show the RMS velocity parameter as a 
function of major axis radius for the five galaxies in our 
sample. We show the inner regions from literature data and 
the outer regions from our data. For all but one galaxy 
(NGC 2768) a general decline in the RMS velocity parameter 
is seen with radius (as might be expected from the generally 
flat rotation velocity and declining velocity dispersion pro- 
files observed). The implications of this for the mass profile 
of each galaxy requires detailed dynamical modeling. 



6 CONCLUSIONS 

We have developed a new technique for extracting spectra of 
the galaxy background contained in the multi-object obser- 
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Figure 20. The Ar parameter (cumulative value of A within 
radius R) with major axis radius scaled by the effective radius. 
The SAURON central values are presented as stars, and our data 
by solid lines. Here NGC 821 is black and NGC 2768 is red. The 
dashed line represents the delimiting value of 0.1 between 'fast-' 
and 'slow- 'rotators (as defined in Emsellem et al. 2007). While 
both galaxies are classified as fast-rotators based on SAURON 
data within 1 R e ff , NGC 821 would be classified as a slow-rotator 
if the data to 3 R e ff are considered. 
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Figure 21. Vrms parameter with major axis radius. Galaxies 
are colour coded, i.e NGC 821 black, NGC 1400 blue, NGC1407 
green, NGC 2768 red and NGC 4494 cyan. The thin lines rep- 
resent the inner literature data, and the thick lines are our data 
out to larger radii (with extrapolated dashed lines connecting the 
two). In general the V rms profiles decline with radius. 

vations using the DEIMOS spectrograph on the Keck tele- 
scope. Here, the slits contained globular clusters, which were 
the primary target of the observing programme. The tech- 
nique allows measurement of the velocity moments (rotation 
velocity, velocity dispersion, and Gauss-Hermite coefficients 
Ii3 and h.4) in a galaxy halo out to at least 3 effective radii 
(R c ff), covering a surface area of some 24 Reff, in compari- 
son to the ~0.4 R^ff achieved by the SAURON integral field 
unit. 

For the 5 early-type galaxies in our sample, we have 
demonstrated that the agreement between our results and 
the literature values is generally extremely good in the re- 
gions of overlap, including the velocity dispersion profiles 
from planetary nebulae data for NGC 821 and NGC 4494. 

For NGC 821 we show that the kinematic signature of 
the disc identified in the surface photometry disappears at 
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around 1.5 R e ff, while the photometry remains highly ellip- 
tical in the outer regions. From this we conclude that NGC 
821 constitutes a highly elliptical bulge with an embedded 
central disc. For NGC 1400, we also detect a falling rota- 
tion profile beyond 1 R e ff. When considered in conjunction 
with the significant I13 found by Spolaor et al. (2008a) in 
the inner regions, this suggests that this early-type galaxy 
also contains an embedded disc. For both NGC 1407 and 
NGC 4494, which are currently classified as pure elliptical 
galaxies, the detection of a trend in the Gauss-Hermite I13 
coefficient suggests that these too harbour a sub-population 
with 'disk-like' kinematics. In NGC 4494, both the 'disk-like' 
kinematics and a coordinated shift in photometric and kine- 
matic position angles provide strong support for the idea 
that this galaxy is a merger remnant. The kinematics of 
NGC 2768 clearly indicates that a disc is present. However, 
in this case, the disc becomes more prominent with galacto- 
centric radius, suggest that this galaxy is, in fact, an SO. It is 
therefore apparent that all five galaxies in our sample show 
evidence for a disc or disc-like components. 

We estimate the angular momentum parameter A_r used 
in the SAURON survey to classify early-type galaxies as ei- 
ther fast or slow rotators. Importantly, our analysis shows 
that the central kinematics are not necessarily representa- 
tive of the kinematics at larger radii, and that galaxies can 
change their rotator class when data from larger radii are 
considered. This has important repercussions for any at- 
tempts to classify galaxies by their central kinematic prop- 
erties alone. 

We have shown that our sample galaxies reveal the rela- 
tively flat rotation and velocity dispersion profiles suggestive 
of massive dark matter halos, thus demonstrating the value 
that our new technique will have when full dynamical mod- 
eling of these galaxies is carried out. 

Finally, we note that the observations used here were 
not optimised for this project. Future work that uses slit 
masks designed for the purpose of extracting galaxy halo 
spectra should provide additional improvements. 
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Parameter 


New value 


Old value 


NGC 1400 


V S y S 


563 ± 5 


560 ± 10 


V rot 


113 ± 6 


139 ± 17 


CO 


281 ± 7 


301 ± 4 


NGC 1407 


V S y S 


1791 ± 5 


1794 ± 10 


V rot 


54 ± 6 


53 ± 19 


(TO 


270 ± 7 


313 ± 9 



Table Al. New and old (Spolaor et al. 2008a) kinematic param- 
eters for NGC 1400 and NGC 1407. 

APPENDIX A: THE RE-ANALYSIS OF 

NGC 1400 AND NGC 1407 LONG-SLIT DATA. 

Spolaor et al. (2008a) used the van der Marel & Franx (1993) 
code to measure velocity moments in 68 and 82 apertures 
along the major axes for NGC 1407 and NGC 1400, re- 
spectively. The spectra had a S/N of >30 A -1 at 5000 A, 
and reached a radial extent of 0.56 R c a (~ 4.11 kpc) for 
NGC 1407 and 1.30 R e s (~ 3.58 kpc) for NGC 1400. Mea- 
surement of moments was then carried out using 17 template 
stars of spectral class F5 to K4. Briefly, the procedure iden- 
tified the template that best fit the central regions (where 
the signal-to-noise is highest) and applied this to all spatial 
slices. The drawback to this method is that the best-fitting 
template was found using only the very central regions and 
the existence of strong stellar population gradients (particu- 
larly metallicity) may require differing templates in the inner 
and outer regions. 

We have therefore re-analysed the data of Spolaor et 
al. (2008a) using the pPXF code of Cappellari et al. (2007). 
The main advantage of this code is that: 

i) A template match is found for all individual spatial 
slices. 

ii) The code finds the combination of templates that 
best matches the data, further improving the fit. 

The results of the new fits are shown in Fig. Al and 
Table Al. Errors on derived parameters were estimated by 
running a series of Monte Carlo simulations. The results do 
not change the general conclusions of Spolaor et al. (2008a). 
However, the kinematically decoupled core identified in Spo- 
laor et al. (2008a) is now much more distinct and the I13 and 
I14 values appear much more consistent. 
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Figure Al. A comparison of the results of the results of Spo- 
laor et al. (2008a) (black points; NGC 1400; top and NGC 1407; 
bottom) to our re-analysis using the pPXF software (red points). 
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